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Abstract 

The short- and intermediate-distance behaviour of the hybrid adiabatic potentials 
is calculated in the framework of the QCD string model. The calculations are per- 
formed with the inclusion of Coulomb force. Spin-dependent force and the so-called 
string correction term are treated as perturbation at the leading potential-type regime. 
Reasonably good agreement with lattice measurements takes place for adiabatic curves 
excited with magnetic components of field strength correlators. 
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§ ! 1 Introduction 

Gluonic degrees of freedom in the nonperturbative region should manifest themselves as QCD 
bound states containing constituent glue, so one expects that purely gluonic hadrons (glue- 
balls) should exist as well as hybrids where the glue is excited in the presence of qq pair. 
General agreement is that the lightest hybrids occur in the mass range between 1.3 and 1.9 
GeV, so the absolute mass scale remains somewhat imprecise in the absence of exact analytic 
methods of nonperturbative QCD. Existing experimental data seem to point towards gluonic 
excitations being present, and prima facie candidates are identified |J, but the conclusive 
evidences have never been presented. There is no hope that in the nearest future data anal- 
yses could unambigiously pin-point the signatures for gluonic mesons and settle the issue of 
constituent glue. The state-of-art is that the predictions of different models on hadronic spec- 
tra and decays are involved in order to distinguish between gluonic mesons and conventional 
ones. 

In such situation the lattice gauge calculations remain the only source of knowledge. 
Lattice calculations are now accurate enough to serve as a guide, so that the results of different 
Q CD-motivated approaches can be compared and contrasted with lattice data. Of particular 
interest are the measurements of gluelump and hybrid adiabatic potentials 0. These 
simulations measure the spectrum of the glue in the presence of static source in the adjoint 
colour representation (gluelump) and in the presence of static quark and antiquark separated 
by some distance R. These systems are the simplest ones and play the role of hydrogen atom 
of soft glue studies, as, first, the gluonic effects are not obscured by light dynamical quarks, 
and, second, the problem of centre-of-mass motion separation is not relevant here. 
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The large distance limit of hybrid adiabatic potentials is important, as one expects the 
formation of confining string at large R. The short range limit is relevant to the heavy hybrid 
mass estimations. One expects that in the case of very heavy quarks the hybrid resides in the 
bottom of the potential well given by the adiabatic curve, which, in accordance with lattice 
results 0, is somewhere around 0.25 fm for lowest curves. 

In the present paper we study hybrid adiabatic potentials in the so-called QCD string 
model This model deals with quarks and point-like gluons propagating in the confining 
QCD vacuum, and is derived from the Vacuum Background Correlators method. In the latter, 
the confining vacuum is parametrised by the set of gauge-invariant field strength correlators 
responsible, among other phenomena, for the area law asymptotics. The basic assumption 
of the QCD string model is the minimal area law for the Wilson loop, so that the only 
nonperturbative input is the string tension a. 

The QCD string model describes the spectra of qq mesons with remarkable agreement 
. |0] . It is also quite successfull in describing glueballs M , hybrids M and gluelump [T(J , as 
well as meson-hybrid-glueball mixing []TT] . 

First studies of hybrid adiabatic potentials in the QCD string model were performed in [|T^] , 
with special attention paid to the large-distance limit. It was shown that at large interquark 
distances two kinds of QCD string vibrations take place, the potential-type longitudinal and 
string-type transverse ones. Here we consider the short-distance behaviour of the excitation 
curves. 

The paper is organised as follows. In Section 2 we briefly discuss the essentials of the 
QCD string approach for gluons. The effective Hamiltonian for a gluon bound by the static 
quark- ant iquark pair is derived in Section 3. It is argued in Section 4 that at short and 
intermediate interquark distances the potential-type regime of string vibrations is adequate, 
and the lowest excitation curves are calculated. The spin-dependent forces and the so-called 
string corrections are considered in Section 5. Results and discussion are given in Section 6 
together with conclusions and outlook. Appendices contain the details of our variational 
calculations. 



2 Gluons in the confining background 

The QCD string model for gluons is derived in the framework of perturbation theory in the 
nonperturbative confining background JIB] . The main idea is to split the gauge field as 

A I1 = B I1 + a M , (1) 

which allows to distinguish clearly between confining field configurations and confined 
valence gluons a^. The valence gluons are treated as perturbation at the confining background. 
We start with the Green function for the gluon propagating in the given external field 

G^x, y) = (D 2 (B)6^ + 2igF, u {B))-\ (2) 
where both covariant derivative D c £ and field strength depend only on the field _B M : 

D c x a (B) = 6 ca d x + gf cba B b x , (3) 
F; V {B) = d,B a v - d v Bl + gr hc B\Bl. (4) 
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Figure 1: Hybrid Wilson loop 



The term, proportional to F^ V {B) in (Q) is responsible for the gluon spin interaction. We 
neglect it for a moment, it will be considered in Section 5. 

Now we use Feynman-Schwinger representation for the quark-antiquark-gluon Green func- 
tion ||, which, for static quark and antiquark, is reduced to the form 

/•OO P 

G(x g ,y g ) = J Q ds J Dz g exp(-K g )(W) B , (5) 

where angular brackets mean averaging over background field. The quantity K g is the kinetic 
energy of gluon (to be specified below), and all the dependence on the vacuum background 
field is contained in the generalized Wilson loop W, depicted in Fig.l, where the contours Tq 
and Tq run over the classical trajectories of static quark and antiquark, and the contour T g 
runs over the gluon trajectory z g in (|). 

The expression (|5|) is the starting point of the QCD string model, as, under the minimal 
area law assumption, the Wilson loop configuration takes the form 

N 2 — 1 

(W) B = exp(-a(S 1 + S 2 )), (6) 

where Si and S 2 are the minimal areas inside the contours formed by quark and gluon and 
antiquark and gluon trajectories correspondingly. 



3 Einbein field form of the gluonic Hamiltonian 

To proceed further we are to fix the gauge in the reparametrization transformations group. 
For the case of static quark and antiquark sources the most natural way to do this is to identify 
the proper time r of the Feynman-Schwinger representation with the laboratory time. Then 
the classical quark and antiquark trajectories are given by 

R. R. 

z Qn = ( r > y) > z Qv = ( r > ~j) > ( 7 ) 

and the action of the system can be immediately obtained from the representation ([5]): 

T ( '2 1 

A = L dT {~2 + ^2~ a So d ^(^ w 'i) 2 -^ w 'i- 

-a £ d(3 2 ^(w 2 w' 2 ) 2 -w 2 2 w'2 2 } ■ (8) 
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Here r is the three-dimensional gluonic coordinate, and the minimal surfaces Si and 5*2 are 
parametrized by the coordinates w^r, /3j), i = 1,2, = -^f L , w' ifM = ^ i£ . Choosing the 
straight-line ansatz for the minimal surfaces one has in the laboratory time gauge 

w i0 = r, wi )2 = ±(l-P)- + /3r. (9) 



The kinetic energy in (|8|) is given in the so-called einbein field form |14| . The einbein field 
/i = /i(r) is the auxiliary field introduced to deal with relativistic kinematics. Note that in 
the case of gluon one is forced to introduce it from the very beginning, as it provides the 
meaningful dynamics for the massless particle. 

In order to pass to the Hamiltonian formulation it is convenient to get rid of Nambu-Goto 
square roots in (§), introducing continuous set of einbein fields z/j = z/j(r, as it was first 
suggested in ||: 

-f^-jf^T' 1 -^ 

'?,2 = <* - 4-(r 1|2 r) 2 , r 1|2 = r±5:. (10 ) 



Note, that the Lagrangian (|H]) describes the constrained system. As no time derivatives 
of the einbeins enter it, the corresponding equations of motion play the role of second-class 



constraints 14 



The Hamiltonian H = pr — L is easily obtained from the Lagrangian fllPl) : 

H = H + ^+ f 1 d{5°-^ + f 1 d(3 2 °-^ + f 1 d(3^ + I' d(3 2 ^, (11) 
2 Jo v\ Jo i>2 Jo 2 Jo 2 



2(/x+ Ji + J 2 ) 

Jl{V + Jl) H n— + J2) 



2A(/x + Ji + J 2 ) I r\ nr lJ r\ 
2J1J2 1 



,2 (rir 2 )(pn)(pr 2 ) (12) 

' 1'2 J 



A = (ji + Jx)(/i + J 2 ) - Ji J 2 ^T-, Ji = /' dPifivitfi), i = l,2. 
The Hamiltonian (O) together with the constraints 



dH 5H 

completely defines the dynamics of the system at the classical level. To quantize one should 
first find the extrema of einbeins from the equations ([H|) and substitute them back to the 
Hamiltonian (p|). Then, the extremal values of einbeins would become the nonlinear operator 
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functions of coordinate and momentum, and, in addition, the problem of operator ordering 
would arise. To avoid this complicated problem the approximate einbein field method is 
usually applied in the QCD string model calculations. Namely, einbeins are treated as c 
number variational parameters: the eigenvalues of the Hamiltonian ([ll]) are found as functions 
of /i and z/j, and minimized with respect to einbeins to obtain the physical spectrum. Such 
procedure, first suggested in ||, provides the accuracy of about 5-10% for the groud state 
(for the details see first entry in j7|). 



4 Potential regime of the QCD string vibrations 



The einbeins /z and Vi(Pi) play the role of constituent gluon mass and energy densities along 
two strings. Note that even with simplifying assumptions of the einbein field method these 
quantities are not introduced as model parameters, but are calculated in the formalism. It is 
clear from the form (|12|) of the kinetic energy that two kinds of motion compete to form the 
spectrum: the potential-type longitudinal with respect to R vibrations due to gluonic mass 
fi and string-type transverse ones due to the string inertia. 

It was shown in [l2|] that for large interquark distances, R ^> l/^/o, these two types 
of motion decouple, displaying the corrections to the leading oR behaviour proportional to 
(■^) 1 / 3 in the case of longitudinal vibrations and proportional to for transverse ones. 

On the contrary, for small R one can neglect the terms Ji responsible for string inertia in 
the kinetic energy (0). Then the Hamiltonian takes the form 



H 



2/i 2 Jo Pl 2z/ x 



2^,2 



or 



d(3 2 — 
o 2u 2 



o 2 



o 2 



(14) 



The spectrum of the Hamiltonian ([14]) was found in |12 |. At small R it reads 



E n (R) = 2 3 / 2 a 1 / 2 (n + 



3A 1 / 2 



o 



3/2 R 2 



23/2( n+ 1)1/2 ' 



(15) 



The extremal values of einbeins are given by 



/i n (i?) = 2 1 /V/ 2 (n + - Al/ " 



( n+ 1)1/20-1/2 



a 



3/2 R 2 



2 l/2 



2 5/2( n+ |)l/2' 

3aV 2 R 2 

2 7 /2( n+ |)V2' 



(16) 
(17) 



where n is the number of oscillator quanta. 

The expressions (|16]), (|T7| ) immediately yield Ji^/n ~ |, so the neglect of string inertia 
is justified. The curves fi n {R) and z/ n (i?) for arbitrary R from |[2f are shown at Fig.|| for 
n = 0,1,2,3 and a = 0.21 GeV. It is clear that the potential-type Hamiltonian can be 
employed at R < 1 fm for n = 0, 1 and at R < 1.5 fm for n = 2, 3, and the corrections due to 
string inertia can be taken into account perturbatively. 

The form ([T4]) allows to eliminate einbeins and arrive at the potential-type Hamiltonian 



H 



+ ari + or 2 . 



(18) 
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Figure 2: Einbein fields v n {R) (solid curve) and fi n {R) (dashed curve) for n = 0,1,2,3 and 
a = 0.21 GeV 2 . 



Nevertheless, as we are going to calculate the spin-dependent forces and string correction, we 
prefer to eliminate only einbeins treating the quantity \i in the framework of the einbein 
field method. 

If only confining force is taken into account, the QCD string model predicts the oscillator 
potential (|I5D with the minimum at R = 0. However, the minimum is shifted, if the long-range 
confining force is augmented by the short-range Coulomb potential, 

V c = + — . 19 

2ri 2r 2 6R K J 



The coefficients in (|1!|) are in accordance with the colour content of the QQg system [ |15 
The QQ Coulomb force in (^) is repulsive, and it is compatible with the behaviour of gluon 
energies || at small R. Note, that such behaviour comes out naturally in the QCD string 
model, as point-like gluon does carry colour quantum numbers. 
The final form of our Hamiltonian reads 

2 

H = f~ + 2 + an + ^ + Vc ' (20) 



The angular momentum is not conserved in the Hamiltonian (p0|) , but it is a good quantum 
number in the einbein field Hamiltonian flH|). For the case of pure confining force we have 
compared the spectra of exact and einbein-field Hamiltonians, and have found that angular 
momentum is conserved within better than 5% accuracy. The same phenomenon is observed 



in the constituent gluon model [fLq| , and seems to be a consequence of linear potential con- 
finement embedded there. 

The eigenvalue problem for the Hamiltonian ([20]) was solved variationally with wave func- 
tions 

^ 7A (r) = Ur) E Ci^Ji^C-)Xx^ (21) 
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Table 1: Quantum numbers of lowest levels 



(a) 


3 PC 


= 1+- 


j = i, l = i, A = 0, 1 


V~ TT 


(b) 


f C 


= l~ 


i = 1, Z = 2, A = 0, 1 




(c) 


f c 


= 2— 


3=2,1 = 2, A = 0,1,2 




(d) 


3 PC 


= 2+- 


3=2, 1 = 3, A = 0,1,2 


TT A 



where ^ is the spin 1 wave function, A = |f | is the projection of total angular momentum 

j onto z axis, z || R. The radial wave functions <f>i(r) were taken to be Gaussian, that is of the 
form exp(— (3 2 r 2 /2) multiplied by the appropriate polynomials, with j3 treated as variational 
parameter. The eigenvalues Eji A {fi,R) = ji A \H (fi, R)\^ji\) were found in such a way, and 
the resulting adiabatic potentials, 



V? lA (R) = E jlA ^*(R),R), 
depend on the extremal value fi* defined from the condition 



dE jlA (n,R) 



0. 



(22) 



(23) 



The details of this variational procedure are given in the Appendix A. 

In the QCD string model the gluon is effectively massive, and has three polarizations 
P|,|10|. Only two of them are excited with magnetic components of field strength correlators, 
used in lattice calculations. We list these states in Table 1, in terms of j,l,A and standard 
notations borrowed from physics of diatomic molecules. For more details justifying such 



correspondence see [10 



Fitting the ground lattice state with Coulomb+linear potential yields the values of pa- 
rameters a, = 0.4 and o = 0.21 GeV 2 . 



5 String correction and spin-dependent interaction 



Now we turn to the calculations of corrections to the leading potential regime (p2|). 

Let us first consider the correction due to string inertia. It corresponds to the terms in 
( JjjD linear in Jj: 

2^ 2 ^ 2^ 2 I r\ 1 + r\ 



= -^V + ^ ^ + ) • (24) 



In the potential regime = ar^ and Jj = OTj/3. So the string correction Hamiltonian takes 
the form 

*- - + ^i) • (*) 

where 

U = r, x p. (26) 

The choice (|26| ) solves the ordering problem in (p5|), as it assures the hermiticity of the 
operator H sc . 
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In actual calculations it is convenient to rewrite (1251) as 



a [1 1\/ , R 2 



where 



HSC -^{{^ + ^{ L+ T H ^ + ^r^ RH ^ (27) 

^ = -(^ + ^ + ^), (28) 

H 2 = zH 1 +pd p d z + d z . (29) 

The adiabatic potentials of string correction are listed in the Appendix B. 

The spin-dependent force originates from the term, proportional to F pu (B) in ([]). This 
term generates the spin- dependent interaction, as 

- iF ik = (SB)«, (30) 

where spin operator S acts at the vector function \l> as 

= -ie ijk V k . (31) 

One averages this term over the background, as it was done in fl7|, and, in our case, it gives 



H mrt = _JL(l_ {LlS) + l_ {Wi) ) t (32) 



HlSm = & [7! {us) + 7I {US) I • (33) 



for nonperturbative and perturbative forces respectively. One easily recognizes the contri- 
bution of Thomas precession in (B^), fl3"3"|). The spin-dependent interaction is conveniently 
represented as 

ffLS<p, ^{fe + 4) (Ls)+ &-4)H' < 35) 

where 

H R = Q -^{-d p + - d^)S + + ^-(d p + - (36) 

jLi P Zj P 

S± = S x i iSy. 
Spin-dependent potentials are given in the Appendix C. 

We would like to stress here, that, in spite of apparently nonrelativistic form of the ex- 
pressions (|2"TD and fl34|), fl3"5|), these are not the nonrelativistic inverse mass expansion. The 
mass p entering these expressions is replaced in matrix elements by the value p*, obtained 
from stationary point equation (|2"3"|). The latter plays the role of effective mass of the gluon, 
and is not large. The -R-dependence of corrections is shown at Figs.pl. f|. 
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Figure 3: Taken with opposite sign potentials of string — V SC (R) (dotted curves), spin- 
dependent perturbative — V LS ^(R) (dashed curves) and nonperturbative — V LS ^ np \R) 
(dashed-dotted curves) corrections along with their sum -V SC (R) - V LS &(R) - V LS ^(R) 
(solid curves), in GeV, for levels (a) and (b) of Table 1. R is measured in fm. 
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Figure 4: The same as in Fig. || for levels (c) and (d). 

6 Results and discussion 

The results of full QCD string calculations are given at Fig.|5| for the excitation curves listed 
in Table 1, together with the lattice data ||. As it is seen from the Fig.|| at small distances 
the calculated curves are in good agreement with lattice results, with the accuracy better 
that 100 MeV. The level ordering is reproduced too, with the exception of S~, Yl g and A g 
levels (Fig.||c). Note, that the lattice data claim only one IT 9 level, and its i?-behaviour is 
rather peculiar (see dashed thick grey curve at Fig.||b). One expects that the curves should 
tend to form degenerate levels as the distance R decreases, fulfilling the angular momentum 
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Figure 5: Hybrid potentials with corrections included (black curves) compared to lattice ones 
(thick grey curves) in units l/r = 400 MeV (r = 2.5 GeV -1 ). QQ-distance R is measured 
in 2r ~ 1 fm. States (a)-(d) are given in Table [l|. Solid, dashed and dashed-dotted curves 
correspond to A = 1,2 and 3. Solid line at right bottom corner of Fig.|5|(a) represents the 
Coulomb+linear potential with a s = 0.4 and a = 0.21 GeV 2 . 



This feature is made explicit in the QCD string model: 
|lT)|l for R = after subtracting the QQ 



conservation demands at R = 0. 
our calculations reproduce the gluelump spectrum 
Coulomb force (last term in ([19|) ) and with obvious replacement 2cr — > o, 



adj 



2a, 



Lattice data indeed seem to follow such tendency. Moreover, the curvatures of all potentials 
but n„ are compatible with the dominance of Coulomb force acting between static sources 



in octet colour representation. Thus one suspects fTSf , fl9|| , that something goes wrong with 



10 



the lattice IT 9 level, and its strange behaviour could be due to the presence of several levels, 
severely mixed and poorly resolved by present simulations. 

The most pronounced feature of the QCD string approach is the following. It was already 
mentioned that the gluon here is effectively massive, and has three polarizations. As a con- 
sequence, the level ordering follows the increasing dimension of valence gluon operator, or, in 
other words, the increasing orbital momentum /. This is in contrast to the standard viewpoint 



of constituent glue studies, see |20| and realization of this idea in the framework of potential 
NRQCD [18]. The level ordering there is supposed to follow the increasing dimension of the 
operators Ei, Bi, DiBk etc. The equations of motion, which relate different operators with 
the same quantum numbers, are involved to exclude spurios states in such picture. 

In our approach, we expect an extra family of levels to appear, namely, and II ff , and 
the corresponding gluelump limit achieved with 1 quantum numbers. The wave functions 
of these states contain mostly the I = component, so these levels should be the lowest 
one-gluon ones. The search for this family, accessible only with electric field correlators, is of 
paramount importance both in gluelump and adiabatic potentials settings. The presence or 
absence of such states would allow to discriminate among models. 



The flux-tube model |2l|, as well as its relativistic version |p2| , assumes that soft glue is 
string-like, with phonon-type effective degrees of freedom. These string phonons are colourless, 
so that the QQ pair is in colour-singlet state. Thus the Coulomb QQ interaction is to be 
attractive. Adiabatic potentials are calculated in f22fl , and, in order to improve the short- 



range behaviour of adiabatic curves, the Coulomb interquark repulsion was added, which 
obviously contradicts the general phylosophy of flux-tube. Without essential modifications 
of dynamical picture at small interquark distances, the flux-tube-type models seem to be 
ruled out by lattice data ||. A rather elaborated constituent gluon model ||I6| , based on 
field-theoretical Hamiltonian approach, agrees with lattice data on hybrid potentials at short 
and intermediate interquark distances only under rather confusing assumption of gluon parity 
taken to be positive. 

The gluelump spectrum, as well as small .R-limit of hybrid adiabatic potentials, is success- 
fully calculated in the bag model p3| . The lowest bag-model gluelump state is l + ~, and the 
X — II splitting at small interquark distances is in accordance with lattice data. In this regard 
we stress once more the importance of lattice measurements with electric field correlators. 
If the ground state £+ — U g family is not found, then, with above-mentioned drawbacks of 
flux-tube and constituent-gluon pictures, it would mean that soft glue is bag-like rather than 
string-like or point-like. 

To conclude, we have presented full QCD string calculations of hybrid adiabatic potentials. 
The results are in general agreement with lattice data. We outline the problems, connected 
with restricted set of gluonic operators used in lattice simulations, and call for further studies 
of excitation curves with electric field operators. 

We are grateful to Yu.A.Simonov for numerous stimulating discussions. The financial 
support of RFFI grants 00-02-17836 and 00-15-96786, and INTAS OPEN 2000-110, 2000-366 
is acknowledged. 

A Variational calculations of adiabatic potentials 

In this Appendix we derive the explicit variational equations for adiabatic potentials 
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The average momentum for oscillatory wave finction with orbital momentum I is given as 

P') r 2 ^, (A.1) 

where (3 is the variational parameter with dimension of mass. Thus the calculation of ex- 
tremum over // (p3[) leads to expression 



S = P\F^- (A.2) 



Due to the symmetry of wave functions, 



Let us introduce the dimensionless variables 



R 

x = (3— fx = f3rx (A.4) 



and define the dimensionless averages 



(fi)=f(x), {±)=g(x). (A.5) 



For the energy levels we will get 



p J K J IV 2 " K ' 12x 

Let us perform now the extremum of the last equation over j3 provided that x is the function 
of and df(x(p))/dp = f'(x)x/(3. So we find that 

fflf \ 2a(f(x)-xf'(x)) 

P (x) = -_= , (A.6) 

J— 3a s (g(x) +xg'(x)) 



E(x) = ^f(x) + P(x) 3a s g(x) + ^ I , (A.7) 

2r 

R(x) = — (A.8) 
P(x) 

The equations ( |A.7| ), ( |A.H| ), along with calculated functions f(x), g(x), P(x) define explicitly 
adiabatic potentials (p2]). Note that all expressions for the averages contain the error function 

erf(x) = 4= I"" dte~ t2 . (A.9) 



SpK JO 
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We present them below for the levels of Table |1|. 
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I) 
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2x e 



7T 
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in 



x I erf(x) ( 1 + — + — + -^=- 

11 4x 2 8x 4 / x/ix 



4x 2 ) 



i Jerf(i) (l + J^j) - - 



-erf (a;) (3 



x + 



2x 
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xe 
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x [ V x 2 x 4 / y 7r x V x 2 
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(fi) 221 = ^ erf(x) 1 + -— - — + - + 



12x 2 8x 4 x 6 / v^rx V3 12x 2 x 



r-4 



1\ 1 [ , / 1 2\ e^ 2 /2x 3 7x 19 4\1 . . « 



2(j f „, x /13 1 6\ e~ x2 /2x 3 llx 7 12" 

x 



< ¥ " + " ^ ~ + ~ + x" + x 3 



2 _ _ _ 

221 — / — 9 

7 3a, f „, s /l 8\ e"* / 35x 2 „ A 8x 6 4x 8 



erf (x) - + — + 16 + — - + 6x 4 + — + — 



2 x 3 1 w V2 x 2 / 0rx I 3 3 3 



4 1 1 \ e"* r 5 I 



= "r f( "H 1 + 3^-4^-4x^J + 7^ ( v 1 + ^ + 2^ 
h/222 x\ KJ \ 2x 2 2x 4 / v/?x V3 x 2 /J v ; 



2(7 (. , /8 1 3 \ e" x / 3 
— erf x )(--_-_ + 4 + — 

x ! V3 x 2 2x 4 / \/tt x \ x 2 



2 

" 2 7 3a s f .. /. 2\ 2 c / . 7x 2 2./- 1 ' 



erf (x) 1 + — - 2 + 



x 3 [ V x 2 / a/ttx \ 3 3 



f , /, 81 33 39 \ e"* /4x 2 1 14 39 \ 1 



2£ I OT .f /„\ ^81 i 66 i 117 \ _ c- x / 117 i 522 i 278a , 68a 3 , 8^ 
7x y ilL \ x ) { 5 -T 5x 2 -T 2x 4 J ^ U 3 5i + 5 ^ 5 ~T 5 

^-if{-6erf(x)(f + ^-^(33 + ^ + ^ + ^ + ^ + ^)} 
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2 23 Ax' 
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I Jp r fMfl 4- - 4- 6 f 22 I M 2 - ^ _ _ 8xM 

^ seri W ^±-h 7a . 4 J -t- 7 ^ ^ x a -t- 15a . 15 15 15 ) 

{«£(*) (f + § " J) " ^ ( * + 
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21 
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7x 3 



{erf(x) (-§ + + ^ (-11 - ^ - ^ + f + ¥ + ^)} 



J / J^9_ 121 13 \ e~ x / 149 13 \ 

(ri)232 ~ T W l 70i 2 140a; 4 + 56xV + J^x \ + 105a; 2 28a; 4 / 
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y'Tl X 
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7/ 232 ~ x 1 61 ^ V ~ 35X 2 + 28a 1 ) + + 15a; + 105/ r 



2a f p , , /129 121 39 \ e- x / 39 177 Vox 
— <^ erf(x) H H =■ H 1 

7x1 U l 5 5a; 2 Ax 4 J yfH V 2a; 3 5a; 15 



232 ~~ ,, ,, / ,.-,/• -,-, ~~ 



6a s \ /36 11 \ e~ x / 106a; 2 52a; 4 8a; 6 * 

71 " 7a^ 1 (X) IT " + 7^ 1^ " "75 15~ " IF, 



B String corrections 

The string correction potentials, Vjf A = (^ji\\H sc \^ji\^, have the following form: 

VSa _ .-{e^^l.il)^^^)}, (B . 1} 



-2 



Vin = ~ 3erf(,)fl + A_iiU^_f5 + f)L (B.2) 



15/3 I v 7 V 2a; 2 6a; 4 / v^Fa; V x 2 
2ax r /28 67 3 



^120 = -T^^erf(a;)(- + — - — ) + 



105/3 r rf ^ y V 3 ' 3a; 2 2a; 4 

e~ x2 ( 3 32 16a; 3 \] 

2ax I c , > ^7 139 3 \ 

e"^ / 3 47 8x 3 \l ,„ . 
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2 " x ( ,, J13 59 3 3 \ 

e ~ X2 f 6 5 17 4x 4a ^\ 
yx 5 2x 3 6x 3 / J ' 

_2<rx J /8 _5_ _3_\ / 3_ J_ _8_^ 

^222 - 21/? jerl(x) + ^ 2 + ^ J + ^[ 2x5 x3 + 3x/ 

_2ax r /81 237 87 9 \ 



i / ■ • - 2 o .,1 ' 



e~ x f 9 18 131 8x 2 8x' 
+ V^x {2X 1 + 35x2 + "35 - + — + "35" J j 

2ax r . /93 277 1 3 \ 
^ = -^A eii(x) ^ + ^-^ + 2x-«) + 
e~ x2 ( 3 8 53 8x 2 8x 4 \l 



v^Fx \ x 4 5x 2 21 2 105/ J ' 

2ax f /129 41 59 3 \ 

'275 I 61 ^ llF + 10x2 ~ 140x4 ~ 8x6/ 



+ 



e~ x2 f 3 47 53 



y^x \4x 4 35x 2 15, 

C Spin-orbit corrections 

Nonperturbative spin-orbit potentials, given as vffi np ^ = ^^aI-H^^I^-m.)) 



bj3x \ V 2x 2 / x 



LS(np) 18(7 f „, s 1 \ e x2 A 8x 2> 



^™ = 35^i erf Wl 1 + ^J-^l 1+ 9. 



,LS(n P ) 2(7 /. 2 \ f „, x /. 3 \ 3e x 



^220^ = I + 3 erf(x) 1- — 



7(3x V x 2 J [ V 2x 2 / a/7tx 



,.2 



^ 221 " 7(3x\ eit{X) {3 2x 2 + x 4 J U + 3x 2 + x^ 
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^ - £{-«(§+si-sO + £(-§ + ?)} 1 (a7) 
vS< " Pl = ik{ ed(x}( - ie+6+ ^~ i ~7r ( - §! + &+4x+ ^}- (a8) 



2e~ x2 / 1 1 16a; 8a; 3 \ 1 _ . 

+-/=" ( 7a TTZ oc 7ac I i ' (C9) 



v^F \x 3 15a: 35 105 

^ (np) = ^{^(M-^+^+^H+m-fk)}- (cio) 



Perturbative spin-dependent potentials, = ( ^jza | -H" ls ^ | ^jiA / 5 read as 



5 [ 2a; 3 ^fix \ x 



7 s \ 5a; 3 Sv^r 



Vm m = -^,P{ -^ + — (6-^ + 4^) \, (CM) 



v-- = 3 ^ (c , 7) 



7 [ x 3 V a: 2 / \3 3a; 2 a; 4 , 

^(p) 8 ferf(x) / 9 X e"*" /18 9 5 12a; 2 8a; 4 \ 1 
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rLS(p) 2 _ ^ferf(x) ( A 3 ^ , e~ x ( 4 3,6 8x r 



v« w = Ti a 'i 3 {-^{- i + ^) + v*{TE--i + ^-l5)\- (a20) 
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